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Abstract. Density is the turbulence statistics that is most readily avail- 
able from observations. Different regimes of turbulence correspond to different 
density spectra. For instance, the viscosity-damped regime of MHD turbulence 
relevant, for instance, to partially ionized gas, can be characterized by shallow 
and very anisotropic spectrum of density. This spectrum can result in sub- 
stantial variations of the column densities. Addressing MHD turbulence in the 
regime when viscosity is not important over the inertial range, we demonstrate 
with numerical simulations that it is possible to reproduce both the observed 
Kolmogorov spectrum of density fluctuations observed in ionized gas by mea- 
suring scintillations and more shallow spectra that are obtained through the 
emission measurements. We show that in supersonic turbulence the high den- 
sity peaks dominate shallow isotropic spectrum, while the small-scale underly- 
ing turbulence that fills most of the volume has the Kolmogorov spectrum and 
demonstrates scale-dependent anisotropy. The limitations of the spectrum in 
studying turbulence induce searches of alternative statistics. We demonstrate 
that a measure called "bispectrum" may be a promising tool. Unlike spectrum, 
the bispectrum preserves the information about wave phases. 



1. Density Fluctuations as a Probe of MHD Turbulence 

The paradigm of interstellar medium has undergone substantial changes recently. 
Instead of quiescent medium with hanging and slowly evolving clouds a turbu- 
lent picture emerged (see review by Ballesteros-Parredes et al. 2006, McKee & 
Ostriker 2007 and ref. therein). With magnetic field being dynamically impor- 
tant and dominating the gas pressure in molecular clouds, this calls for studies 
of compressible magnetohydrodynamic (MHD) turbulence. 

Key ideas in describing MHD turbulence can be traced back to Iroshnikov 
(1963) and Kraichnan (1965) work and the classical work that followed (see 
Montgometry & Turner 1981, Higdon 1984, Montgomery, Brown & Matthaeus 
1987). A more recent model by Goldreich & Sridhar (1995, henceforth GS95) has 
been successfully tested by numerical simulations (Cho & Vishniac 2000, Maron 
& Goldreich 2001, Cho, Lazarian & Vishniac 2002, Cho & Lazarian 2002, 2003, 
henceforth CL03). 

Density statistics is the easiest to infer from observations. In comparison, 
velocity spectra require elaborate techniques to be used (see Lazarian 2006 and 
ref. therein). How informative can be density fluctuations for understanding 
of MHD turbulence is the subject of the present review. However, we should 
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mention that for many important astrophysical applications, e.g. for interstel- 
lar chemistry, for star formation, for propagation of radiation etc. the density 
fluctuations themselves are crucially important. 

Subsonic compressible MHD is rather well studied topic today. It is sugges- 
tive that there may be an analogy between the subsonic MHD turbulence and its 
incompressible counterpart, namely, GS95 model. Therefore the correspondence 
between the the two revealed in CL03 is expected. 

It could be easily seen, that in the low-beta case density is perturbed mainly 
due to the slow mode (CL03). Slow modes are sheared by Alfven turbulence, 
therefore they exhibit Kolmogorov scaling and GS95 anisotropy for low Mach 
numbers. However, for high Mach numbers we expect shocks to develop. Density 
will be perturbed mainly by those shocks. 

One can also approach the problem from the point of view of underlying 
hydrodynamic equations. It is well known that there is a multiplicative sym- 
metry with respect to density in the ideal flow equations for an isothermal fluid 
(see e.g. Passot & Vazquez-Semadeni, 1998). This assume that if there is some 
stochastic process disturbing the density it should be a multiplicative process 
with respect to density, rather than additive, and the distribution for density 
values should be lognormal, rather that normal. The aforementioned work shows 
that for ID numerical simulations of high-Mach hydrodynamics the distribution 
is approximately lognormal, having power-law tails in case of 7 7^ 1. 

In MHD the above described symmetry is broken by the magnetic tension. 
However, numerical studies in Beresnyak, Lazarian <fe Cho (2005, henceforth 
BLC05) show an approximate correspondence of the PDFs obtained with MHD 
simulations to the lognormal scaling. 

With a high sonic Mach we expect a considerable amount of shocks arise. In 
a sub- Alfvenic case, however, we expect oblique shocks be disrupted by Alfvenic 
shearing, and, as most of the shocks are generated randomly by driving, almost 
all of them will be sheared to smaller shocks. The evolution of the weak shocks 
will be again governed by the sonic speed, and structures from shearing as in 
low Mach case should arise. 

We also note that shearing will not affect probability density function (PDF) 
of the density, but have to affect its spectra and structure function (SF) scaling. 
In other words, we deal with two distinct physical processes, one of which, 
random multiplication or division of density in presence of shocks, affect PDF, 
while the other, Alfvenic shearing has to affect anisotropy and scaling of the 
structure function of the density. The numerical studies that we describe below 
confirm the theoretical considerations above. 



2. Viscosity-Damped MHD Turbulence as a Special Case 

Turbulence can be viewed as a cascade of energy from a large injection energy 
scale to dissipation at a smaller scale. The latter is being established by equating 
the rate of turbulent energy transfer to the rate of energy damping arising, for 
instance, from viscosity. Naively, one does not expect to see any turbulent 
phenomena below such a scale. 

Such reasoning may not be true in the presence of magnetic field, however. 
Consider magnetized fluid with viscosity v much larger than magnetic diffusivity 
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r/, which is the case of a high magnetic Prantl number Pr fluid. The partially 
ionized gas can serve as an example of such a fluid up to the scales of ion-neutral 
decoupling (see a more rigorous treatment in Lazarian, Vishniac & Cho 2004, 
henceforth LVC04). Fully ionized plasma is a more controversial example. For 
instance, It is well known that for plasma the diffusivities are different along and 
perpendicular to magnetic field lines. Therefore, the plasma the Prandtl number 
is huge if we use the parallel diffusivity un S> u±. A treatment of the fully ionized 
plasma as a high Prandtl number medium is advocated in Schechochihin et al. 
(2004, henceforth SCTMM). 

The turbulent cascade in a fluid with isotropic rj proceeds up to a scale at 
which the cascading rate, which for the Kolmogorov turbulence, i.e. vi ~ I 1 / 3 , 
is determined by the eddy turnover rate v\/l gets equal to the damping rate 
r]/l 2 . Assuming that the energy is injected at the scale L and the injection 
velocity is Vl, the damping scale l c is Li?e _3//4 , where Re is the Reynolds num- 
ber LVl/t]. The corresponding scale l c varies from ~ 10 -2 pc for Warm Neutral 
Medium to 10 -4 for molecular clouds (LVC04). It is evident, however, that mag- 
netic fields at the scale l c at which hydrodynamic cascade would stop are still 
sheered by eddies at the larger scales. This should result in creating magnetic 
structures at scales <C l c (LVC04). Note, that magnetic field in this regime is 
not a passive scalar, but important dynamically. Thus, GS95 turbulence range 
an inertial range in the partially ionized gas from the injection scale to l c , but 
below l c is the domain of the viscosity-damped turbulence, where plasma/gas 
compressed within current sheets can contribute to the observed extreme scatter- 
ing events and the formation of the Small Ionized or Neutral Structures (SINS) 
(Lazarian 2007). Most compressions arises in the situations when the magnetic 
fields exhibits perfect reversals. 

Cho, Lazarian & Vishniac (2002) reported a new viscosity-damped regime 
of turbulence using incompressible MHD simulations (see an example of more 
recent simulations in Fig. 1). In this regime, unlike hydro turbulence, motions, 
indeed, do not stop at the viscosity damping scale, but magnetic fluctuations 
protrude to smaller scales. Interestingly enough, these magnetic fluctuations 
drive small amplitude velocity fluctuations at scales <C l c - CL03 confirmed 
these results with compressible simulations and speculated that these small 
scale magnetic fluctuations can compress ambient gas to produce ubiquitous 
tiny structures observed in ISM (see also Fig. lb). According to the model in 
LVC04, while the spectrum of averaged over volume magnetic fluctuations scales 
as Es(k) ~ k^ 1 , the pressure within intermittent magnetic structures increases 
with the decrease of the scale as (Sb)^ ~ k, while the filling factor 4>k ~ k , the 
latter being consistent with numerical simulations. 

In Fig. [1] we show the results of our high resolution compressible MHD 
simulations that exhibit strong density fluctuations at the scales below the one 
at which hydrodynamic turbulence would be damped. According to LVC04 
model, the viscosity-damped regime is ubiquitous in turbulent partially ionized 
gas. Some of its discussed consequences, such as an intermittent resumption of 
the turbulence the fluid of ions as magnetic fluctuations reach the ion-neutral 
decoupling scale, are important for radio scintillations. 
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Figure 1. Left: Filaments of density created by magnetic compression in 
the slice of data cube of the viscosity-damped regime of MHD turbulence. 
Right: Spectra of density and magnetic field are similar, while velocity is 
damped. The resistive scale in this regime is not L/Rm but LRmT 1/2 .(tr om 
Beresnyak & Lazarian, in preparation) . 



3. Simulations 3D statistics and Simulated Observations 
3.1. Spectrum 

The power spectrum (PS) of density fluctuations is an important property of a 
compressible flow. In some cases, the spectrum of density can be derived analyt- 
ically. For nearly turbulent motions in the presence of a strong magnetic field, 
the spectrum of density scales similarly to the pressure, i.e. E p {k) ~ k- 7 / 3 if 
we consider the polytropic equation of state p = ap 1 (jH). In weakly magnetized 
nearly incompressible MHD turbulence, however, velocities convect density fluc- 
tuations passively inducing the spectrum E p (k) ~ /c -5 / 3 |22|). In supersonic 
flows, these relations are not valid anymore because of shocks accumulating 
matter into the local and highly dense structures. Due to the high contrast of 
density, the linear relation 5p = c^Sp is no longer valid, and the spectrum of den- 
sity cannot be related to pressure so straightforwardly. In addition, the strong 
asymmetry of density fluctuations suggests the need to analyze the logarithm of 
density instead of density itself. 

In Figure [2] we present the PS of fluctuations of density for models with 
different M s . As expected, we note a strong growth of the amplitude of density 
fluctuations with the sonic Mach number at all scales. This behavior is observed 
both in sub-Alfvenic as well as in super- Alfvenic turbulence (see Fig. [2]). In 
Table [T] we calculate the spectral index of density a p and the logarithm of density 
aiogp within the inertial ranges estimated from the PS of velocity. The width 
of the inertial range is shown by the range of solid lines with slopes —5/3 and 
—2 in all spectra plots. It is estimated to be within k £ (3,20). In Table [T] we 
also show the errors of estimation which combine the error of the fitting of the 
spectral index at each time snapshot and the standard deviation of variance of 
«p,iogp in time. The slopes of the density spectra do not change significantly 
with M s for subsonic experiments and correspond to analytical estimations 
(about —2.2, which is slightly less than —7/3, for turbulence with Ma ~ 0.7 and 
about —1.7, which is slightly more than —5/3, for weakly magnetized turbulence 
with Ma ~ 7.3). Such an agreement confirms the validity of the theoretical 
approximations. Those, nevertheless, do not cover the entire parameter space. 
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Figure 2. Spectra of density for experiments with different values of A4 S 
and with Ma ~ 0.7 (left panel) and Ma ~ 7.3 (right panel) for models with 
medium resolution (256 3 ). Grey error bars signify the variance of spectra with 
time. The solid lines with slopes —5/3 and —2 cover the inertial range used 
to estimate the spectral indices of p spectra (from UH). 



anisotropy for logarithm of density onisotropy for logarithm of density 




Figure 3. Left: Anisotropy for the 2™ d -order SF for the logarithm of density 
for Ma ~ 7.3; Right: the same is for Ma ~ 7.3. The logarithm of density 
traces the anisotropies of turbulent magnetic field. 



While the fluid motions become supersonic, they strongly influence the density 
structure, making the small-scale structures more pronounced, which implies 
flattening of the spectra of density fluctuations (see values for A4 S > 1.0 in 
Table [T]) (see also BLC05). Although, the spectral indices are clearly different 
for sub- and super-Alfvenic turbulence, their errors are relatively large (up to 
±0.3). These errors have been calculated by summing the maximum value of the 
uncertainty of fits for individual spectra at each time snapshot and the standard 
deviation of time variation of spectral indices. The uncertainty of fit contributes 
the most to the total error bar (about 60%- 70%). 



3.2. Anisotropies 

Another question is the anisotropy of density and the logarithm of density 
structures. For subsonic turbulence it is natural to assume that the density 
anisotropies will mimic velocity anisotropies in GS95 picture. This was con- 
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Table 1. Slopes of the Power Spectrum of Density and the Logarithm of 
Density Fluctuations (from [12) 
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The values have been estimated within the inertial range for models with Ma ~ 0.7 (left) 
and Ma ~ 7 (right). Errors of spectral indices combine the errors of estimation at each time 
snapshot and the standard deviation of variance in time. Errors for sonic Mach numbers are 
the standard deviation of their variance in time calculated over the period starting from t > 5 

to the last available snapshot. 



firmed in CL03, who, however, observed that for supersonic turbulence the con- 
tours of density isocorrelation get round, corresponding to isotropy. BLC05, 
however, showed that anisotropies restore the GS95 form if instead of density 
one studies the logarithm of density. This is due to the suppression of the influ- 
ence of the high density peaks, which arise from shocks. It is these peaks that 
mask the anisotropy of weaker, but more widely spread density fluctuations. 

In Figure ?? we show lines that mark the corresponding separation lengths 
for the second-order SFs parallel and perpendicular to the local mean magnetic 
fieldQ In the case of subAlfvenic turbulence, the degree of anisotropy for density 
is very difficult to estimate due to the high dispersion of points. However, rough 
estimates suggest more isotropic density structures, because the points extend 
along the line h\ ~ For models with Ma ~ 7.3 the points in Figure ?? 
have lower dispersion, and the anisotropy is more like the type from GS95, i.e. 

2/3 , 

Zjj ~ l_l . In both Alfvenic regimes, the anisotropy of density does not change 
significantly with A4 S . Plots for the logarithm of density show more smooth 
relations between parallel and perpendicular SFs. The dispersion of points is 
very small. Moreover, we note the change of anisotropy with the scale. Lower 
values of SFs correspond to lower values of the separation length (small-scale 
structures), so we might note that the logarithm of density structures are more 
isotropic than the GS95 model at small scales, but the anisotropy grows a bit 
larger than the GS95 prediction at larger scales. This difference is somewhat 
larger in the case of models with stronger external magnetic field (compare plots 
in the left and right columns of Figure ??), which may signify their dependence 
on the strength of B ext . The anisotropy of log p structures is marginally de- 
pendent on the sonic Mach number, similar to the density structures. All these 



1 The local mean magnetic field was computed using the procedure of smoothing by a 3D Gaus- 
sian profile with the width equal to the separation length. Because the volume of smoothing 
grows with the separation length I, the direction of the local mean magnetic field might change 
with I at an arbitrary point. This is an extension of the procedures employed in @. 
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Figure 4. The bispectra for density 
ki and k2- Here we compare different 
turbulence. The structure is clearly 
supersonic (right panel) cases. 
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show the degree of correlation between 
sonic numbers for strongly magnetized 
different for subsonic (left panel) and 



observations allow us to confirm the previous studies (see BLC05) that suggest 
that the anisotropy depends not only on the scale but on Ma- 

3.3. Bispectrum 

Attempts to use multipoint statistics are a more traditional way to remove the 
constraints that the use of two point statistics, e.g. power spectra entails. Un- 
fortunately, very high quality data is needed to obtain the multipoint statistics. 
Among multipoint statistics, bispectrum (see Scoccimarro 1997) seems the most 
promising. This is partially because it has been successfully used in the studies 
of the Large Scale Structure of the Universe. 

Bispectrum is a Fourier transform of the three point correlation function 
and if the power spectrum P(k) is defined as 

<Mki)<Wk 2 )> = ^(k)fe(k! + k 2 ) (1) 

where 5d is the Dirac delta function that is zero apart from the case when 
ki + k 2 = 0, the bispectrum B123 is 

(Spik^S p(k 2 )Sp(k 3 )) = B 123 5 D (ki + k 2 + k 3 ). (2) 

The bispectra of density shown in Fig. [5] give information as to how shocks 
and magnetic fields effect turbulence. It has been shown by (fl2l ) and (@) that in 
supersonic turbulence shocks produce compressed density and a shallower spec- 
trum. Looking at Fig. 0] it is clear that these shocks play a crucial role in the 
correlation of modes. The subsonic cases show little correlation between any 
points except the case of k\ = k%. The compressed densities from supersonic 
turbulence are critical for correlations between frequencies since waves get closer 
together and therefore have a much higher interaction rate. We find an good 
agreement with (fl2l ) in that the density structures are affected by the presence 
of magnetic field. When there is a weak magnetic field present it is clear from 
Figure |4] that the system lacks the stronger correlations that are characteris- 
tic of the sub-Alfvenic models. For super-Alfvenic simulations correlations in 
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frequency are not as readily made due to large dispersion of density structure. 
The bispectrum of supersonic hydro models are similar to the super-Alfvenic, 
supersonic cases. A measure related to bispectrum is bicohence (see Koronovskii 
& Hramov 2002). 

4. Summary 

Above we discussed the statistics of density in compressible MHD turbulence. 
We analyzed spectra, structure functions for experiments with different sonic 
A4 S and Alfven Ma Mach numbers. Our results are as follows: 

• The viscosity-damped regime of MHD turbulence relevant to partially 
ionized gas can be characterized by shallow and very anisotropic spectrum of 
density. This spectrum can result in large variations of the column densities. 

• The amplitude of density fluctuations strongly depends on M s both in 
weakly and strongly magnetized turbulent plasmas. 

• The flattening observed in density spectra is due to the contribution of 
the highly dense small-scale structures generated in the supersonic turbulence. 

• Fluctuations of the logarithm of density are much more regular than those 
of density. The logarithm of density exhibits the GS95 scalings and anisotropies. 

• Bispectra show strong correlations for supersonic models than for subsonic 
ones. It suggests the importance of shocks in mode correlations. 

Acknowledgement of the support of the NSF Center for Magnetic Self-Organization in 
Laboratory and Astrophysical Plasma. 
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